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Abstract 

Time-series  measurements  of  velocity,  temperature  and  conductivity  on  the  northern  California 
shelf  during  two  winter  seasons  permit  an  observational  test,  in  vertically  integrated  form,  of  a 
simple  set  of  subinertial  momentum  and  heat  balances  for  the  bottom  boundary  layer,  which  have 
resulted  from  recent  theoretical  work.  These  are:  (1)  an  along-isobath  momentum  equation  that 
reduces  to  a  classical  Ekman  balance;  (2)  a  cross-isobath  momentum  equation  in  which  the  Ekman 
balance  is  modified  by  a  buoyancy  force  caused  by  distortion  of  the  isopycnal  surfaces  within 
the  boundary  layer;  and  (3)  a  heat  balance  in  which  variability  of  temperature  is  produced  by 
cross-isobath  advection.  The  measurements  confirm  the  importance  of  buoyancy  in  the  cross¬ 
isobath  momentum  equation,  and,  as  has  recently  been  predicted  theoretically,  they  indicate  that 
buoyancy  is  a  dominant  effect  when  the  boundary  layer  is  thick,  which  typically  occurs  during 
downwelling-favorable  flows.  An  Ekman  balance  describes  subinertial  fluctuations  in  the  along- 
isobath  momentum  equation  with  moderate  success.  The  mean  along-isobath  momentum  equation 
is  significantly  influenced  by  a  buoyancy  force  caused  by  an  along-isobath  temperature  gradient, 
and  along-isobath  advection  is  as  important  as  cross-isobath  advection  in  the  heat  balance.  Thus 
along-isobath  variability  of  temperature,  which  has  been  neglected  in  most  theoretical  work,  is 
fundamental  to  the  structure  of  the  velocity  field  and  the  evolution  of  the  temperature  field. 

The  most  striking  feature  of  the  measurements  is  a  persistent  cross-isobath  mean  flow,  which 
is  offshore  below  20  m  above  bottom,  onshore  above  20  m,  and  uniform  above  30  m.  The  transport 
produced  by  this  flow,  relative  to  the  unsheared  flow  above  30  m,  is  much  larger  than  the  Ekman 
transport  associated  with  the  mean  along-isobath  bottom  stress.  Thus  the  dominant  force  pro¬ 
ducing  a  mean  cross— isobath  transport,  relative  to  the  unsheared  flow  above  30  m,  is  not  bottom 
stress,  but  must  instead  be  some  other  mechanism.  A  likely  explanation  is  a  vertically  nonuniform 
along-isobath  density  gradient,  which  is  poorly  resolved  by  the  measurements. 
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1  Introduction 

Recent  theoretical  work  on  the  bottom  boundary  layer  has  focused  on  a  problem  in  which  a  subin- 
ertial,  along-isobath  flow  begins  from  rest  in  a  stably  stratified  ocean  above  a  gently  sloping,  non¬ 
conducting  sea  floor  (see  the  review  by  Garrett  et  al.  [1993]  and  the  recent  articles  by  Ramsden 
[1995a,  1995b]  and  Middleton  and  Ramsden  [1996]).  In  this  problem,  the  sloping  bottom  permits 
the  cross-isobath  Ekman  transport  to  modify  the  density  field,  and  the  combination  of  the  sloping 
bottom  and  the  distortion  of  isopycnal  surfaces  by  mixing  and  advection  produces  cross-isobath 
buoyancy  forces  that  influence  the  velocity  field.  The  interaction  between  density  and  velocity  leads 
to  evolution  and  structure  that  are  profoundly  different  from  the  more  classical  results  for  a  flat 
bottom.  In  particular,  if  the  forcing  is  steady,  the  models  predict  evolution  toward  a  steady  state 
in  which  the  along-isobath  velocity  is  brought  to  zero  at  the  bottom  by  a  frictionless  thermal-wind 
balance  that  occurs  within  the  boundary  layer,  so  that  the  bottom  stress  and  cross— isobath  Ekman 
transport  vanish  and  the  bottom  offers  no  frictional  resistance  to  the  overlying  flow. 

Although  idealized,  the  recent  theoretical  work  provides  a  concise  set  of  vertically  integrated, 
subinertia!  momentum  and  heat  balances  that  can,  at  least  in  principle,  be  tested  observationally. 
These  are:  (1)  an  along-isobath  Ekman  balance,  in  which  bottom  stress  is  proportional  to  cross¬ 
isobath  transport;  (2)  a  cross-isobath  momentum  equation  in  which  the  Ekman  balance  is  modified 
by  a  buoyancy  force;  and  (3)  a  heat  balance  in  which  cross— isobath  advection  produces  temporal 
variability  of  temperature.  These  balances  represent  a  simple  but  oceanographically  relevant  set  of 
boundary  layer  processes,  and  they  are  fundamental  to  the  existing  understanding  of  the  bottom 
boundary  layer  and  its  response  to  the  overlying  flow.  In  particular,  the  concept  of  an  along-isobath 
Ekman  balance  is  a  cornerstone  of  geophysical  fluid  dynamics  (e.g.,  Pedlosky,  1979). 

The  above  balances  have  not,  to  our  knowledge,  been  tested  observationally.  In  fact,  few  oceanic 
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studies  have  resolved  the  structure  of  the  bottom  boundary  layer.  Several  researchers  have  measured 
velocity  and  temperature  within  a  few  meters  of  the  bottom,  without  corresponding  measurements 
in  the  upper  part  of  the  boundary  layer  (see,  for  example,  the  reviews  by  Grant  and  Madsen  [1986], 
Dyer  and  Soulsby  [1988]  and  Cacchione  and  Drake  [1990]).  Other  workers  have  reported  coarsely 
spaced  measurements  in  the  upper  part  of  the  boundary  layer,  without  corresponding  measure¬ 
ments  near  the  bottom  (Kundu,  1976;  D’Asaro,  1982a;  D’Asaro,  1982b;  Lentz  and  Trowbridge, 
1991).  Observations  with  profiling  instrumentation  have  resolved  the  entire  boundary  layer,  but 
the  durations  have  been  too  short  for  a  definitive  examination  of  the  sub-inertial  dynamics  (Mer¬ 
cado  and  van  Leer,  1976;  Weatherly  and  Martin,  1978;  Armi  and  D’Asaro,  1980;  Dickey  and  van 
Leer,  1984).  Recent  observations  on  continental  slopes  have  focused  on  the  structure  up  to  hundreds 
of  meters  above  bottom,  and  they  have  not  resolved  the  much  thinner  frictional  boundary  layer  in 
which  Ekman  dynamics  are  relevant  (Thorpe,  1987;  Thorpe  et  al.,  1990;  White,  1994).  Weatherly’s 
(1972)  moored  array  resolved  the  bottom  boundary  layer  in  the  Florida  current,  but  the  duration 
(one  week)  was  too  short  for  a  statistically  significant  examination  of  subinertial  variability,  and 
the  analysis  was  limited  primarily  to  estimating  bottom  stress,  determining  the  thickness  of  the 
logarithmic  layer,  and  documenting  Ekman  veering.  Bird  et  al.  (1982)  reported  measurements  that 
spanned  the  bottom  boundary  layer  on  the  Bermuda  Rise  for  a  period  of  eight  months,  but  the 
vertical  resolution  was  coarse  (four  sensors  between  0.8  and  62  m  above  bottom)  and  the  analysis 
was  limited  primarily  to  a  qualitative  comparison  with  simulations  based  on  a  turbulence  closure 
model. 

We  report  measurements  obtained  as  part  of  the  Sediment  Transport  Events  on  Shelves  and 
Slopes  (STRESS)  program,  which  resolve  the  vertical  structure  of  the  bottom  boundary  layer  on  the 
northern  California  shelf  during  each  of  two  winter  seasons.  The  purpose  of  the  analysis  is  to  test 
the  vertically  integrated  momentum  and  heat  balances  that  have  resulted  from  recent  theoretical 
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work.  In  contrast  to  previous  observations,  the  STRESS  measurements  have  a  duration  sufficiently 
long  to  examine  subinertial  variability  and  they  span  the  boundary  layer  with  vertical  resolution 
sufficiently  fine  to  estimate  the  terms  in  the  vertically  integrated  balances.  In  the  following,  we 
first  present  the  theoretical  background  (Section  2)  and  describe  the  measurements  and  analysis 
(Section  3).  We  then  present  and  discuss  the  results  (Sections  4  and  5)  and  summarize  conclusions 
(Section  6). 

2  Background 

2.1  Mathematical  model 

The  mathematical  model  addresses  an  idealized  configuration  in  which  a  turbulent  boundary  layer 
above  a  non-conducting,  gently  sloping  sea  floor  is  overlain  by  a  stratified  interior  in  which  turbulent 
fluxes  of  momentum  and  heat  are  negligible  (Figure  1).  The  coordinate  system  is  defined  so  that 
the  x-y  plane  coincides  with  the  sloping  sea  floor,  with  x  cross-isobath,  y  along-isobath,  and  z 
perpendicular  to  the  bottom,  with  z  —  0  at  the  bottom.  Time  is  denoted  t.  The  time  scales  of 
interest  and  the  along-isobath  and  cross-isobath  spatial  scales  are  assumed  to  be  sufficiently  large 
that  local  and  advective  accelerations  are  negligible  in  comparison  with  the  Coriolis  acceleration. 
The  relationship  between  density  p  and  temperature  T  is  assumed  to  be  approximately  linear; 
i.e.,  p  ~  po  —  P(T  —  To),  where  po  is  a  fixed  reference  density  and  /?  and  To  are  constants.  The 
temperature  gradients  dT/dx  and  dT/dy  are  assumed  to  be  approximately  independent  of  z  within 
the  boundary  layer.  This  configuration  is  essentially  the  same  as  in  the  theoretical  work  reviewed 
by  Garrett  et  al.  (1993),  although  it  is  slightly  more  general  in  that  dT/dy  is  not  assumed  to  be 
zero. 

The  dynamics  of  the  boundary  layer  are  described  by  simplified  momentum  and  heat  equations 
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for  a  Boussinesq  fluid  (e.g.,  Pedlosky,  1979).  The  x  and  y  momentum  equations  reduce  approx¬ 
imately  to  a  balance  between  Coriolis  acceleration,  pressure  gradient  and  stress  divergence,  with 
an  additional  gravitational  term  resulting  from  the  sloping  coordinate  system.  The  x  momentum 
equation  is 

-  ^ = -to + %  ■  ■  (i) 


and  the  y  momentum  equation  is 


pofu  = 


dp  dljy 
dy  dz 


The  z  momentum  equation  is  approximately  hydrostatic: 


(2) 


0  = 


(3) 


and  the  heat  equation  reduces  approximately  to  a  balance  between  time  dependence,  cross-isobath 
and  along-isobath  advection,  and  vertical  mixing: 


dT  dT  dT  d<j> 
dt  U  dx  V  dy  dz' 


(4) 


Here  g  is  gravitational  acceleration,  /  is  the  Coriolis  parameter,  a  is  the  small  bottom  slope,  (u,  v,  w) 


is  the  velocity  vector,  p  is  pressure,  rzx  and  rzy  are  components  of  the  stress  tensor,  <f>  represents 
turbulent  flux  of  temperature,  and  T,  p  and  p0  are  respectively  temperature,  density  and  a  fixed 
reference  density,  as  before. 

In  (4),  the  term  w  dT/dz  has  been  neglected  without  obvious  justification.  In  the  theoretical 
work  reviewed  by  Garrett  et  al.  (1993),  this  term  is  identically  zero  because,  by  assumption,  the 
velocity  is  every where  parallel  to  the  bottom,  so  that  w  =  0.  Here  we  point  out  that  even  if  w 
is  not  precisely  zero,  w  dT/dz  is  small  in  comparison  with  udT/dx  within  the  boundary  layer, 
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provided  that  the  temperature  in  the  boundary  layer  is  relatively  well  mixed.  To  see  this,  note 
that  w  is  at  most  O(au),  except  in  front-like  regions  with  short  cross-isobath  or  along-isobath 
scales,  which  are  excluded  from  the  analysis.  In  addition,  dT/dx  is  0{adToo/dz ),  where  is  the 

interior  temperature  (e.g.,  Garrett  et  al.,  1993).  Thus  u  dT/dx  is  0(audT0O/dz )  while  wdT/dz  is 
0(au  dT/dz ),  so  that  w  dT/dz  <  u  dT/dx  provided  that  dT/dz  <  dT^/dz.  Above  the  boundary 
layer,  where  the  temperature  is  not  well  mixed,  neglect  of  w  dT /dz  is  not  necessarily  justified. 

Neglect  of  temporal  accelerations  in  (1)  and  (2)  is  a  conspicuous  omission,  particularly  in  (2), 
because  shelf  measurements  typically  indicate  that  dv/dt  is  comparable  to  fu  for  subinertial  fluctu¬ 
ations  (e.g.,  Allen,  1980).  Neglect  of  temporal  accelerations  is,  however,  an  acceptable  approxima¬ 
tion  in  the  present  analysis.  Because  a  proper  treatment  of  temporal  acceleration  is  complicated, 
the  justification  is  postponed  until  the  discussion  (Section  5.3). 


2.2  Vertically  integrated  momentum  balances 


To  obtain  a  set  of  vertically  integrated  balances,  it  is  convenient  to  define  z  =  6  to  be  a  position  just 
above  the  turbulent  boundary  layer,  where  turbulent  fluxes  of  heat  and  momentum  are  negligible. 
The  quantity  <5  varies  with  time  as  the  turbulent  boundary  layer  thickens  and  thins. 


Begin  by  differentiating  (1)  and  (2)  with  respect  to  z,  substituting  the  linear  relationship  be¬ 
tween  p  and  T,  and  using  (3)  to  eliminate  p.  The  result  is  a  pair  of  equations  describing  the  2 


derivatives  of  u  and  v: 


-09 


(dT  dT\ 

U*  adz) 


d2rzx 
+  dz2  ’ 


(5) 


and 


d2Tzy 
dz 2  ' 


(6) 


Note  that  the  terms  involving  temperature  in  (5)  and  (6)  represent  the  familiar  effect  of  thermal 
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wind,  expressed  in  a  coordinate  system  in  which  the  x  axis  is  inclined  at  a  small  angle  with  respect 
to  horizontal. 

Next,  integrate  (5)  and  (6)  from  an  arbitrary  value  of  z  to  z  =  <5,  where  drzx/dz  and  drzy/dz 
are  zero.  The  resulting  expressions  are 

-/»»/(» =  +  (7) 

and 

p0f(u  -us)  =  (3g  jf  ^  dz'  +  ,  (8) 

where  subscript  S  denotes  evaluation  at  z  =  6  and  z1  is  a  dummy  variable  of  integration.  Equations 
(7)  and  (8)  are  a  restatement  of  (1)  and  (2),  with  quantities  involving  p  and  p  expressed  in  terms 
of  u$,  vs  and  T. 

By  integrating  (7)  and  (8)  from  z  =  0  to  z  —  6,  one  obtains  a  pair  of  vertically  integrated 
momentum  balances: 


rd  6 2  dT  rd 

-p0fjo(v-vs)dz  =  l3g  -—  +  ajo(T-Ts)dz 


7~bx  > 


(9) 


and 


f0  62  dT 

Pa S  jo  (“  -  “()  dz  =  ~  Th 


(10) 


where  Tbx  and  r^y  are  the  x  and  y  components  of  the  bottom  stress.  The  first  and  last  terms  in 
(9)  and  (10)  represent  the  classical  Ekman  balance.  The  terms  involving  temperature  represent 
buoyancy  forces  that  result  from  non-horizontal  isopycnal  surfaces.  In  the  special  case  in  which 
dT/dy  is  zero,  (9)  and  (10)  are  consistent  with  the  dynamics  in  the  theoretical  work  reviewed  by 
Garrett  et  al.  (1993),  and  they  are  identical  to  the  expressions  used  in  the  vertically  integrated 
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model  developed  by  Trowbridge  and  Lentz  (1991). 


2.3  Vertically  integrated  heat  balance 

A  simple  form  of  the  vertically  integrated  heat  balance  can  be  obtained  by  integrating  (4)  from 
z  =  0  to  z  =  6  and  applying  the  condition  $  —  0  at  the  non-conducting  sea  floor  and  in  the 
frictionless  interior.  The  result  is 


fs  8T  ,  dT  fs  ,  dT  fs  , 

l  miz=-^L%iz-^Lviz- 


(ii) 


The  term  on  the  left  side  of  (11)  is  the  vertically  integrated  temporal  variability  of  temperature 
within  the  boundary  layer.  The  terms  on  the  right  side  represent  the  effects  of  cross-isobath  and 
along-isobath  advection,  respectively.  Like  (9)  and  (10),  (11)  is  consistent  with  the  theoretical 
work  reviewed  by  Garrett  et  al.  (1993)  in  the  special  case  in  which  dT/dy  is  zero. 

A  more  dynamically  relevant  form  of  the  heat  balance  involves  the  temperature  anomaly  T  — 
because  this  quantity  appears  in  the  buoyancy  force  in  (9)  and  because  it  is  related  to  stratification, 
which  influences  vertical  mixing.  Because  a  heat  balance  involving  T  -  Ts  is  substantially  more 
complicated  than  (11)  and  requires  special  treatment,  we  postpone  consideration  of  this  balance 
until  Section  5. 


3  Methods 

3.1  Measurements 

The  STRESS  observations  (Fredericks  et  al.,  1993;  Trowbridge  and  Nowell,  1994)  occurred  on 
the  northern  California  shelf  (Figure  2)  in  approximately  the  same  location  as  the  Coastal  Ocean 
Dynamics  Experiment  (CODE;  see  Beardsley  and  Lentz,  1987).  The  shelf  in  the  STRESS  region 
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consists  of  a  narrow,  steeply  sloping  inner  portion,  between  the  coast  and  the  70-m  isobath,  and 
a  broader,  gently  sloping  outer  portion,  between  the  70-m  isobath  and  the  shelf  break,  at  a  depth 
of  approximately  150  m.  The  primary  site,  denoted  C3,  was  at  a  depth  of  approximately  90  m.  A 
secondary  site  on  the  same  isobath,  denoted  C3',  was  displaced  from  C3  by  approximately  5  km 
to  the  southeast.  An  additional  secondary  site,  denoted  C4,  was  at  a  depth  of  approximately  130 
m,  approximately  8  km  offshore  of  C3.  We  use  a  coordinate  system  in  which  the  y  axis  is  positive 
toward  a  bearing  of  317°  relative  to  true  north,  which  coincides  with  the  isobath  orientation  at 
C3,  and  the  x  axis  is  positive  toward  a  bearing  of  47°.  The  topography  near  C3  has  well  defined 
cross-isobath  and  along-isobath  directions,  and  our  results  axe  not  sensitive  to  small  changes  in  the 
orientation  of  the  coordinate  system. 

During  1988-89  (STRESS-1),  A.  J.  Williams,  of  the  Woods  Hole  Oceanographic  Institution, 
deployed  a  bottom  tripod  at  C3,  and  B.  Butman,  of  the  U.  S.  Geological  Survey,  deployed  subsurface 
moorings  at  C3  and  C3'  (Figure  3).  The  tripod  supported  six  BASS  acoustic  current  meters 
and  eight  temperature  sensors  between  heights  of  0.2  and  6  m  above  bottom.  The  subsurface 
mooring  at  C3  supported  five  vector  averaging  current  meters  (VACMs),  which  record  velocity 
and  temperature,  at  heights  between  6  and  30  m,  and  the  subsurface  mooring  at  C3'  supported 
VACMs  at  6  and  18  m  above  bottom.  Three  of  the  VACMs  at  C3  were  fitted  with  temperature- 
conductivity  sensors  and  optical  transmissometers,  and  both  VACMs  at  C3'  were  fitted  with  optical 
transmissometers.  The  bottom  tripod  was  within  300  m  of  the  subsurface  mooring  at  C3.  The 
tripod  and  moorings  were  on  site  from  6  December  1988  to  24  January  1989  (STRESS-la)  and 
from  26  January  1989  to  27  February  1989  (STRESS-lb).  All  of  the  instrumentation  functioned 
satisfactorily  during  both  of  these  periods. 

STRESS-1  coincided  with  the  Surface  Mixed  Layer  Experiment  (SMILE;  see  Alessi  et  al.  [1991] 
and  Dever  and  Lentz  [1994]).  As  part  of  SMILE,  S.  J.  Lentz  and  R.  C.  Beardsley,  of  the  Woods 
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Hole  Oceanographic  Institution,  deployed  a  surface  mooring  at  C3.  This  mooring  supported  vector 
measuring  current  meters  (VMCMs),  which  record  velocity  and  temperature,  between  a  depth  of 
50  m  and  the  water  surface.  We  use  the  SMILE  measurements  briefly  (Section  4.1)  to  illustrate 
the  qualitative  structure  above  the  STRESS  array. 

During  1990-91  (STRESS-2),  A.  J.  Williams  deployed  bottom  tripods  and  B.  Butman  deployed 
subsurface  moorings  at  both  C3  and  C4  (Figure  3).  At  each  site,  the  tripod  was  within  600  m 
of  the  mooring.  Each  tripod  supported  BASS  current  meters  and  temperature  sensors  between 
heights  of  0.1  and  6  m  above  bottom.  The  mooring  at  C3  supported  an  array  of  VACMs  fitted 
with  temperature-conductivity  sensors  and  transmissometers,  as  well  as  an  array  of  temperature- 
conductivity  sensors.  The  mooring  at  C4  supported  a  VACM,  an  array  of  VACMs  fitted  with 
temperature-conductivity  sensors  and  transmissometers,  a  temperature-conductivity  sensor,  and 
an  array  of  temperature  sensors.  Both  tripods  and  both  moorings  were  on  site  and  successfully 
collecting  data  from  12  January  1991  to  8  March  1991.  Two  of  the  BASS  sensors  at  each  site,  one  of 
the  VACM  velocity  sensors  at  C3,  and  one  of  the  VACM  conductivity  sensors  at  C3  malfunctioned 
(Figure  3). 

All  of  the  STRESS  instruments  obtained  measurements  essentially  continuously  but  recorded 
at  a  number  of  different  rates.  In  the  analysis,  we  use  measurements  that  have  been  averaged  and 
resampled  to  produce  time  series  of  hourly  averaged  values. 

3.2  Estimation  of  the  terms  in  the  vertically  integrated  balances 

The  STRESS  measurements  provide  hourly  averaged  time  series  of  u,  v  and  T  with  relatively 
fine  vertical  resolution  at  C3  and  coarser  resolution  at  C3'  and  C4.  To  estimate  the  terms  in 
the  vertically  integrated  balances,  we  use  the  measurements  at  C3  to  evaluate  quantities  involving 
vertical  structure,  and  we  use  the  measurements  at  C3/  and  C4  to  provide  information  about  cross- 
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isobath  and  along-isobath  gradients.  To  focus  on  sub-inertial  dynamics,  we  first  estimate  the  terms 
in  the  vertically  integrated  balances  by  using  the  hourly  averaged  data,  and  we  then  apply  filter  pl64 
(e.g.,  Limeburner,  1985),  which  leaves  the  energetic  subinertial  fluctuations  nearly  unattenuated, 
while  effectively  removing  the  energetic  tidal  and  near-inertial  fluctuations  (Section  4.1). 

Some  aspects  of  this  procedure  are  straightforward.  Appropriate  values  of  p0,  f  and  a  are 
1000  kg  m~3,  9.1  x  10-5  s_1,  and  0.005,  respectively  (Figure  2).  To  determine  /?  we  use  the 
best-fit  values  obtained  from  the  temperature-conductivity  measurements,  which  are  0.35  kg  m-3 
deg-1  C  for  STRESS-1  and  0.29  kg  m-3  deg-1  C  for  STRESS-2  (Figure  4).  To  put  the  velocity 
and  temperature  measurements  at  C3  on  a  common  grid  suitable  for  evaulating  vertical  integrals, 
we  use  the  heights  of  the  functioning  current  sensors  as  nodes.  For  nodes  at  which  temperature 
measurements  were  not  obtained,  we  assign  a  temperature  record  by  interpolating  between  the 
nearest  adjacent  temperature  measurements.  This  operation  results  in  a  grid  of  eleven  nodes  for 
STRESS-1  and  eight  nodes  for  STRESS-2.  We  use  a  midpoint  rule  to  estimate  integrals  with 
respect  to  z  and  centered  differences  to  estimate  time  derivatives. 

To  determine  S,  we  assume  that  the  turbulent  boundary  layer  coincides  with  the  bottom  mixed 
layer.  We  therefore  set  6  equal  to  the  height  of  the  lowest  sensor  above  the  mixed  layer,  except 
in  the  relatively  rare  occasions  when  the  mixed  layer  extends  to  the  top  of  the  STRESS  array, 
when  we  set  S  equal  to  the  height  of  the  top— most  sensor.  We  estimate  mixed— layer  thickness 
by  using  the  procedure  described  by  Lentz  and  Trowbridge  (1991),  in  which  the  mixed  layer  is 
defined  to  extend  to  the  highest  sensor  at  which  the  local  temperature  is  within  AT  of  the  bottom¬ 
most  temperature,  where  AT  is  an  estimate  of  the  uncertainty  in  the  temperature  measurements. 
Here  we  use  AT  =  0.05°  C,  as  in  Lentz  and  Trowbridge  (1991).  Use  of  mixed-layer  thickness  to 
estimate  the  thickness  of  the  turbulent  boundary  layer  is  debatable,  because  some  measurements 
of  turbulence  energy  dissipation  indicate  that  the  two  layers  do  not  always  coincide  (Dewey  et  al., 
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1988).  Use  of  mixed-layer  thickness  is  simple,  however,  and  it  is  necessary  in  our  case  because 
we  have  no  measurements  of  turbulence  statistics.  For  the  most  part,  estimates  of  mixed-layer 
height  are  below  the  height  of  the  top-most  sensor  (Section  4.2,  Figures  12  and  13),  so  that  the 
measurements  resolve  the  entire  mixed  layer.  Exceptions  occur  primarily  during  the  first  half  of 
STRESS-2  (Section  4.2,  Figure  13),  when  the  mixed  layer  reached  the  top  sensor  during  a  number 
of  events,  so  that  the  top  of  the  mixed  layer  may  in  fact  have  been  above  the  top  of  the  measurement 
array. 

Direct  estimation  of  dT/dx  is  possible  for  STRESS-2  but  not  for  STRESS-1.  To  obtain  a  time 
series  of  dT/dx  during  STRESS-2,  we  use  the  mean  of  the  estimates  at  nine  heights  between  0.19 
and  19  m  based  on  temperature  differences  between  C3  and  C4.  The  resulting  time  series  of  dT/dx 
has  a  mean  of  0.087°  C  per  km  and  a  standard  deviation  of  0.045°  C  per  km.  The  time  series  of 
dT/dx  at  the  different  heights  are  remarkably  consistent  (temporal  means  and  standard  deviations 
vary  by  less  than  10%  at  the  different  heights  and  the  squared  correlation  coefficient  r 2  is  greater 
than  0.67  for  all  pairs  of  estimates).  This  consistency  justifies  the  assumption  (Section  2)  that 
dT/dx  is  independent  of  z.  To  estimate  dT/dx  during  STRESS-1,  we  use  a  fixed  value  equal  to 
the  mean  of  the  time  series  obtained  during  STRESS-2.  Although  crude,  this  procedure  captures 
the  dominant  feature  of  the  bottom  temperature  field,  which  is  the  offshore  cooling  associated  with 
increasing  depth.  The  STRESS-1  results  are  insensitive  to  the  precise  value  of  dT/dx  that  we  use, 
and  the  STRESS-2  results  are  changed  only  slightly  if  we  use  a  fixed  value  of  dT/dx  rather  than 
the  measured  time  series. 

Direct  estimation  of  dT/dy  is  possible  for  STRESS-1  but  not  for  STRESS-2.  To  determine 
dT/dy  during  STRESS-1,  we  use  the  mean  of  the  estimates  at  two  heights  (6  and  18  m)  based  on 
temperature  differences  between  C3  and  C37.  The  resulting  time  series  of  dT/dy  has  a  mean  of 
-0.015  °  C  per  km  and  a  standard  deviation  of  0.020°  C  per  km.  The  time  series  of  dT/dy  at  the  two 
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heights  are  only  moderately  well  correlated  (r2  =  0.37),  but  they  have  the  same  dominant  features, 
which  are  a  roughly  constant  negative  value  with  brief  periods  of  positive  values.  To  estimate 
dT/dy  during  STRESS-2,  we  use  a  fixed  value  equal  to  the  mean  of  the  time  series  obtained 
during  STRESS-1.  This  approach  is  probably  inaccurate,  since  the  standard  deviation  of  dT/dy 
during  STRESS-1  is  larger  than  the  mean,  and  since  the  STRESS-1  results  for  the  heat  balance 
change  significantly  if  we  use  a  fixed  value  of  dT/dy  instead  of  the  measured  time  series.  However, 
use  of  a  fixed,  representative  value  at  least  permits  us  to  evaluate  the  importance  of  dT/dy  during 
STRESS-2. 

To  estimate  bottom  stress,  we  use  the  log-profile  method  (e.g.,  Sternberg,  1968),  which  is  based 
on  the  standard  unstratified  wall-layer  model  (e.g.,  Monin  and  Yaglom,  1971),  in  which  the  flow 
speed  varies  with  z  according  to 


\/v?  +  v2  =  —  ln(^)  -  —  ln(^o) ,  (12) 

AC  AC 

and  the  flow  coincides  in  direction  with  the  bottom  stress.  Here  is  the  shear  velocity,  k  is  von 
Karman’s  empirical  constant,  approximately  equal  to  0.40,  and  zq  is  the  bottom  roughness  length. 
To  apply  the  method,  one  fits  velocity  measurements  to  (12),  using  u*/k  and  («*/ k)  ln(zo)  as  free 
parameters,  and  then  determines  the  magnitude  of  the  bottom  stress  from  the  definition  |  rj  |=  pou2. 
Here  we  fit  velocity  measurements  obtained  by  the  bottom  three  functioning  BASS  sensors  to  (12), 
and  we  assume  that  the  bottom  stress  coincides  in  direction  with  the  velocity  averaged  over  the 
bottom  three  sensors.  We  use  three  sensors  because  the  logarithmic  model  applies  only  at  small 
heights  and  because  three  is  the  minimum  number  that  permits  an  assessment  of  the  accuracy 
of  the  stress  estimates  based  on  discrepancies  in  the  model  fits.  In  addition,  the  lowest  complex 
empirical  orthogonal  function  (EOF;  see,  for  example  Davis  [1976];  here  based  on  the  complex 
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velocity  u  +  iv),  which  accounts  for  more  than  85%  of  the  variance  in  the  hourly  velocity  data  from 
the  entire  current  meter  array,  indicates  that  while  the  speed  profiles  are  logarithmic  throughout 
the  BASS  array,  the  direction  is  vertically  uniform  (particularly  during  STRESS- 1)  only  at  the 
lowest  three  sensors,  with  subtle  but  consistent  changes  in  direction  at  the  fourth  sensor  and  above. 
Thus  the  model  (12),  which  implies  a  logarithmic  speed  profile  and  a  vertically  uniform  direction, 
appears  to  apply  only  at  the  bottom  three  sensors.  The  relationship  between  shear  velocity  and 
flow  speed  that  results  from  the  log-profile  estimates  is  roughly  consistent  with  a  quadratic  drag 
law  based  on  zo  =  0.01  m  (Figure  5),  which  corresponds  to  a  drag  coefficient  q  =|  t%  |  /[po(u2  +  v2)] 
of  7.5  x  10~3  at  z  =  1  m. 

Log-profile  estimates  of  bottom  stress  are  uncertain  because  of  uncertainties  in  current  mea¬ 
surements  and  because  of  shortcomings  in  the  logarithmic  model  caused,  for  example,  by  distortions 
of  the  near-bottom  velocity  produced  by  acceleration  (Soulsby  and  Dyer,  1981)  and  stratification 
due  to  temperature  or  suspended  sediment  (e.g.,  Monin  and  Yaglom,  1971).  To  quantify  the  effect 
of  errors  in  current  measurements  in  a  simple  way,  assume  for  the  moment  that  the  model  (12) 
is  correct,  that  errors  in  current  measurements  have  a  Gaussian  probability  distribution,  and  that 
the  standard  deviation  of  measurement  errors  is  equal  to  the  standard  deviation  of  the  ensemble  of 
all  residuals  from  fits  of  (12)  to  the  hourly  averaged  velocity  measurements  from  the  bottom  three 
BASS  sensors,  which  is  0.0014  m/s  for  STRESS-1  and  0.0030  m/s  for  STRESS-2.  Given  these 
assumptions,  the  95%  confidence  intervals  for  estimates  of  are  ±0.001  m/s  for  STRESS-1  and 
±0.002  m/s  for  STRESS-2,  indicating  uncertainties  in  stress  estimates  that  are  much  smaller  than 
the  stresses  themselves  (Figure  5).  Our  opinion  is  that  uncertainties  in  «*  estimates  are  in  fact  sub¬ 
stantially  larger  than  ±0.001  or  ±0.002  m/s,  and  that  they  are  caused  primarily  by  shortcomings  in 
the  logarithmic  model.  However,  calculations  of  the  effect  of  flow  acceleration  based  on  the  model 
of  Soulsby  and  Dyer  (1981),  the  effect  of  stratification  due  to  temperature  based  on  an  approach 
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similar  to  that  recently  proposed  by  Friedrichs  and  Wright  (1997),  and  the  effect  of  stratification 
due  to  suspended  sediment  based  on  the  model  developed  by  Glenn  and  Grant  (1987)  indicate  that 
these  effects  were,  for  the  most  part,  important  only  when  the  bottom  stress  was  weak,  so  that  they 
are  not  important  for  the  present  purposes.  More  importantly,  with  a  single  exception  (Section 
5.1),  failures  to  close  the  vertically  integrated  momentum  balances  did  not  occur  during  periods  in 
which  any  of  these  effects  appeared  to  be  large. 

4  Results 

4.1  Qualitative  temporal  variability  and  vertical  structure 

Before  proceeding  to  the  vertically  integrated  balances,  a  brief  examination  of  temporal  variability 
and  vertical  structure  provides  a  useful  qualitative  characterization  of  the  observations.  Here  we 
present  rotary  spectra,  filtered  time  series,  vertical  profiles  of  statistics,  and  estimates  of  gradient 
Richardson  number. 

Rotary  spectra  indicate  that  velocity  fluctuations  occurred  predominantly  in  a  subinertial  band 
(absolute  value  of  frequency  less  than  approximately  0.5  cpd),  a  diurnal  band,  a  near— inertial  band, 
and  a  semidiurnal  band  (Figure  6).  Diurnal  and  semidiurnal  fluctuations  had  predominantly  clock¬ 
wise  polarization,  corresponding  to  negative  frequencies  in  Figure  6,  while  subinertial  currents  were 
primarily  along-isobath  and  therefore  had  no  particular  polarization.  During  STRESS-1,  approxi¬ 
mately  55%  of  the  velocity  variance  at  the  top-most  sensor  was  subinertial,  and  approximately  35% 
occurred  in  the  combination  of  diurnal,  semidiurnal  and  near-inertial  bands.  In  contrast,  approxi¬ 
mately  85%  of  the  variance  at  the  top-most  sensor  during  STRESS-2  was  subinertial.  During  both 
STRESS-1  and  STRESS-2,  fluctuations  in  bottom  stress  were  predominantly  subinertial  (Figure 
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Time  series  of  filtered  (i.e.,  subinertial)  along-isobath  velocity  indicate  energetic  fluctuations 
with  time  scales  of  several  days,  during  which  the  flow  was  predominantly  poleward  with  brief 
periods  of  equatorward  flow  (Figures  8a  and  9a).  Time  series  of  filtered  temperature  (Figures 
8b  and  9b)  indicate  long  periods  with  strong  stratification  and  other  periods  when  temperature 
throughout  the  measurement  array  was  well  mixed.  The  qualitative  variability  of  temperature  is 
generally  consistent  with  a  simple  model  based  on  the  first  two  terms  in  (11),  in  which  poleward 
flow,  corresponding  to  downwelling  Ekman  transport,  produces  warming,  and  equatorward  flow, 
corresponding  to  up  welling  Ekman  transport,  produces  rapid  cooling.  However,  there  are  obvious 
exceptions  to  this  simple  view  (e.g.,  day  40  during  STRESS-1).  Thick  mixed  layers  generally 
correspond  to  periods  of  strong  poleward  flow,  as  reported  by  Lentz  and  Trowbridge  (1991). 

Statistics  of  velocity  and  temperature  provide  a  useful  means  of  illustrating  vertical  structure. 
The  mean  cross-isobath  velocity  for  both  STRESS-1  and  STRESS-2  is  offshore  near  the  bottom 
and  onshore  above  approximately  20  m  above  bottom,  with  a  uniform  structure  at  heights  above 
30  m  (Figure  10a).  The  magnitude  and  structure  of  the  cross-isobath  mean  flow  are  remarkably 
consistent,  especially  considering  the  small  magnitude  and  the  fact  that  three  different  sensors 
(BASS,  VACM  and  VMCM)  produced  the  measurements.  The  mean  along-isobath  flow  is  poleward 
at  all  sensors,  with  a  magnitude  generally  much  larger  than  that  of  the  cross-isobath  mean  flow 
(Figure  10b).  The  standard  deviations  of  filtered  along-isobath  and  cross-isobath  velocities  are 
comparable  to  the  means  (i.e.,  0.01  to  0.03  m/s  for  u  and  0.10  to  0.15  m/ s  for  u).  The  lowest  complex 
EOF  (based  on  the  complex  velocity  u  +  iv),  which  accounts  for  more  than  93%  of  the  variance 
in  the  filtered  velocity  fields,  indicates  a  vertical  structure  in  which  the  magnitude  (Figure  10c)  is 
strongly  sheared  near  the  bottom  and  the  orientation  (Figure  lOd)  is  consistent  in  both  magnitude 
and  sense  with  simple  models  of  the  unstratified  turbulent  Ekman  layer  (e.g.,  Ellison,  1956).  The 
mean  temperature  profiles  indicate  similar  stable  stratification  during  STRESS— 1  and  STRESS— 
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2,  although  temperatures  were  larger  during  STRESS-2  (Figure  lOe).  Fluctuations  in  filtered 
temperature  have  a  standard  deviation  of  roughly  0.5°  C,  and  the  lowest  EOF  of  filtered  temperature 
fluctuations,  which  accounts  for  more  than  95%  of  the  variance  in  the  filtered  temperature  fields, 
indicates  relatively  large  fluctuations  near  the  bottom  during  STRESS-1  and  a  nearly  uniform 
vertical  structure  during  STRESS-2  (Figure  lOf). 

Calculations  of  gradient  Richardson  number  Ri,  based  on  the  hourly  velocity  and  temperature 
measurements  and  defined  by  Ri  =  (3g(dT/dz)/[(du/dzf  +  (dv/dz)%  provide  a  qualitative  indica¬ 
tion  of  the  likely  importance  of  stable  stratification  in  suppressing  vertical  mixing.  Histograms  of  Ri 
based  on  finite  difference  calculations  for  both  STRESS-1  and  STRESS-2  indicate  values  generally 
well  below  the  critical  value  of  1/4  indicated  by  linear  theory  for  onset  of  instability  (Miles,  1961; 
Howard,  1961)  within  the  bottom  two  meters,  and  values  generally  greater  than  1/4,  with  rare 
excursions  below  1/4,  at  heights  above  roughly  ten  meters  (Figure  11).  Between  heights  of  one  and 
ten  meters,  the  modal  value  of  Ri  gradually  increases,  and  above  fifteen  meters  the  distribution  of 
Ri  is  approximately  independent  of  height.  These  calculations  support  the  use  of  the  unstratified 
wall  layer  model  (12)  to  estimate  bottom  stress  based  on  measurements  within  roughly  two  meters 
of  the  bottom.  The  calculations  of  Ri  indicate  that  stable  stratification  probably  had  a  significant 
effect  on  vertical  mixing  at  heights  greater  than  a  few  meters  above  bottom. 

4.2  Vertically  integrated  momentum  balances 

To  examine  the  vertically  integrated  momentum  balances,  we  regard  the  terms  on  the  right  sides 
of  (9)  and  (10)  as  forcing,  and  we  regard  the  terms  on  the  left  sides  as  response.  One  measure 
of  the  performance  of  the  vertically  integrated  momentum  balances  is  the  squared  correlation 
coefficient  between  forcing  and  response,  denoted  r2.  A  standard  analysis  based  on  Gaussian 
statistics  indicates  that  for  STRESS-1  the  critical  value  of  r2  required  to  reject  the  null  hypothesis 
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that  r2  =  0  at  the  95%  confidence  level  is  approximately  0.21  for  (9)  and  0.13  for  (10).  The 
corresponding  critical  values  of  r2  for  STRESS-2  are  0.33  and  0.16,  respectively.  For  both  STRESS- 
1  and  STRESS-2,  these  estimates  correspond  to  a  temporal  correlation  scale  (i.e.,  the  ratio  of  record 
length  to  effective  number  of  degrees  of  freedom)  of  roughly  4  days  for  (9)  and  2  days  for  (10). 

We  first  consider  the  cross-isobath  momentum  balance.  For  both  STRESS-1  and  STRESS-2, 
the  mean  and  standard  deviation  of  the  response,  —pof  fo(v~v5)  dz>  are  larger  than  the  mean  and 
standard  deviation  of  the  forcing  term  in  a  simple  Ekman  balance,  —  Ux,  by  a  factor  of  roughly 
two  (Table  1).  In  addition,  —pof  fo(v  —  v$)dz  is  not  significantly  correlated  with  — at  the 
95%  confidence  level  (Figure  12a  and  13a).  Thus  the  measurements  do  not  support  a  conventional 
Ekman  balance  in  the  cross-isobath  momentum  equation.  In  contrast,  the  standard  deviation 
of  the  buoyancy  force  in  (9)  is  comparable  to  the  standard  deviation  of  -pof  fo(v  -  v$)dz,  and 
the  mean  of  the  buoyancy  term  is  sufficient,  together  with  the  bottom  stress,  to  close  the  mean 
cross-isobath  momentum  balance  to  within  35%  (Table  1),  which  is  remarkably  good  agreement, 
considering  the  magnitude  of  the  fluctuations.  In  addition,  incorporation  of  the  buoyancy  term 
improves  the  correlation  between  forcing  and  response  to  r2  ~  0.55  (Figures  12b  and  13b),  which 
is  easily  significant  at  the  95%  confidence  level.  Much  of  the  improvement  in  r 2  is  probably  due 
simply  to  the  fact  that  both  the  response  and  the  buoyancy  force  tend  to  increase  in  magnitude  as 
6  increases.  However,  the  fact  that  the  forcing  captures  the  mean  response  and  the  magnitude  of 
the  fluctuations  in  response  indicates  that  the  good  agreement  in  Figures  12b  and  13b  is  physically 
meaningful,  rather  than  being  an  artifact  of  dependence  on  6.  Thus  the  measurements  indicate 
clearly  that  buoyancy  is  an  important  part  of  the  cross-isobath  momentum  balance. 

For  both  STRESS-1  and  STRESS-2,  the  buoyancy  force  in  the  cross-isobath  momentum  equa¬ 
tion  is  large  when  the  mixed  layer  is  thick  (Figures  12  and  13).  In  STRESS-1,  the  buoyancy  force  is 
much  larger  than  the  along-isobath  bottom  stress  during  two  brief  downwelling  events  (i.e.,  events 
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with  offshore  bottom  stress)  centered  on  days  -12  and  40  (Figure  13b).  During  STRESS-2,  the 
buoyancy  force  is  much  larger  than  the  bottom  stress  during  two  downwelling  events  of  longer 
duration,  centered  on  days  25  and  35,  and  also  during  a  brief  upwelling  event  (i.e.,  an  event  with 
onshore  bottom  stress)  centered  on  day  50.  For  the  most  part,  these  observations  are  consistent 
with  theoretical  work  (e.g.,  Trowbridge  and  Lentz,  1991)  indicating  that  mixed  layers  are  thick 
during  downwelling-favorable  flows,  and  that  a  thick  mixed  layer  produces  a  large  cross-isobath 
buoyancy  force  that  becomes  the  dominant  effect  balancing  -pof  fo(v  ~  vs)  dz. 

We  next  consider  the  along-isobath  momentum  balance.  In  this  case,  the  standard  deviation 
of  the  response,  p0f  Jq(u  -  us)dz,  is  comparable  to  the  forcing  term  in  a  simple  Ekman  balance, 
— Tby,  and  both  are  much  larger  than  the  standard  deviation  of  the  buoyancy  term  (Table  1).  The 
correlation  between  forcing  and  response  (r2  ~  0.33)  is  easily  significant  at  the  95%  confidence 
level  (Figures  14a  and  15a),  although  it  is  smaller  than  the  corresponding  correlation  for  the  cross¬ 
isobath  momentum  balance.  Addition  of  the  buoyancy  term  produces  a  marginal  and  probably 
statistically  insignificant  improvement  in  the  correlation  between  forcing  and  response  only  for 
STRESS-2  (Figures  14b  and  15b).  The  buoyancy  term  is  important  in  the  mean  along-isobath 
momentum  balance,  since  the  mean  buoyancy  force  is  approximately  half  of  the  mean  along-isobath 
bottom  stress  (Table  1).  The  mean  along-isobath  momentum  balance  closes  within  20%  for  both 
STRESS-1  and  STRESS-2  (Table  1).  The  presence  of  a  dynamically  significant  buoyancy  force 
in  the  along-isobath  momentum  balance  is  not  consistent  with  the  theoretical  work  reviewed  by 
Garrett  et  al.  (1993),  in  which  along-isobath  variability  was  assumed  to  be  zero. 

It  is  noteworthy  that  the  comparisons  of  forcing  and  response  in  the  momentum  balances  are 
dominated  by  a  few  large  events,  particularly  in  the  case  of  the  cross-isobath  momentum  equation 
(Figures  12  and  13).  Thus  the  effective  number  of  degrees  of  freedom  is  small,  in  spite  of  the  fact 
that  the  record  length  is  much  longer  than  in  previous  observational  studies  that  resolve  the  vertical 
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structure  of  the  bottom  boundary  layer. 

4.3  Vertically  integrated  heat  balance 

To  examine  the  vertically  integrated  heat  balance,  we  regard  the  term  on  the  left  side  of  (11)  as 
response  and  the  terms  on  the  right  side  as  forcing.  As  in  the  case  of  the  momentum  balances, 
one  measure  of  the  performance  of  the  vertically  integrated  heat  balance  is  the  squared  correlation 
coefficient  r2  between  forcing  and  response.  A  standard  analysis  based  on  Gaussian  statistics 
indicates  that  for  (11)  the  critical  value  of  r2  required  to  reject  the  null  hypothesis  that  r2  =  0 
at  the  95%  confidence  level  is  0.06  for  STRESS-1  and  0.07  for  STRESS-2.  Both  of  these  values 
correspond  to  a  temporal  correlation  scale  (ratio  of  record  length  to  effective  number  of  degrees  of 
freedom)  of  roughly  1  day. 

Although  statistically  significant  at  the  95%  confidence  level,  the  correlation  between  response 
and  forcing  for  the  heat  balance  (11)  is  small  if  we  consider  only  cross-isobath  advection  (Figures  16a 
and  17a).  The  correlation  between  forcing  and  response  improves  dramatically  if  we  include  along- 
isobath  advection  (Figures  16b  and  17b).  The  standard  deviations  of  cross-isobath  and  along- 
isobath  advection  are  nearly  identical,  and  the  mean  along-isobath  advection  is  substantially  larger 
than  the  mean  cross-isobath  advection  (Table  1).  As  in  the  momentum  balance,  the  importance 
of  along-isobath  advection  in  the  heat  balance  is  not  consistent  with  the  theoretical  work  reviewed 
by  Garrett  et  al.  (1993)  in  which  along-isobath  variability  was  neglected. 

Note  that  the  mean  heat  balance  does  not  close  well,  particularly  for  STRESS-2  (Table  1).  For 
both  STRESS-1  and  STRESS-2,  the  estimates  of  the  advective  terms  in  the  heat  balance  indicate 
that,  in  the  mean,  the  boundary  layer  should  be  warming  more  rapidly  than  is  observed.  This 
discrepancy  is  likely  a  result  of  the  fact  that  the  measurements  do  not  resolve  dT jdx  and  dT / dy 
particularly  well. 
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5  Discussion 

5.1  Discrepancies  in  the  momentum  balances 

In  STRESS-1,  discrepancies  in  the  cross-isobath  momentum  balance  occurred  primarily  during 
brief  events  near  day  -10  and  just  before  day  40  (Figure  12)  and  discrepancies  in  the  along-isobath 
momentum  balance  occurred  during  the  same  periods  and  also  near  day  10  (Figure  14).  An  exam¬ 
ination  of  these  periods  is  instructive. 

During  day  -10,  transmissometer  records  indicate  the  largest  turbidities  that  were  recorded 
during  STRESS-1,  apparently  because  of  wave-induced  sediment  resuspension  during  a  severe 
storm.  Stable  stratification  by  suspended  sediment  might  have  influenced  the  near-bottom  flow 
during  this  period,  resulting  in  over-estimation  of  the  magnitude  of  bottom  stress  by  the  log-profile 
method.  This  idea  is  consistent  with  analyses  by  Gross  et  al.  (1992)  and  Wiberg  et  al.  (1994) 
of  optical  transmission  measurements  obtained  during  the  STRESS  program,  which  suggest  that 
sediment  concentrations  were  large  enough  during  a  few  severe  storms  to  influence  the  near-bottom 
flow.  Also  in  support  of  this  idea,  estimates  of  stress  obtained  from  a  quadratic  drag  law  with  a 
fixed  drag  coefficient  based  on  z0  =  0.01  m  remove,  essentially  completely,  the  discrepancies  during 
this  period  in  both  (9)  and  (10).  Use  of  a  quadratic  drag  law  with  a  fixed  drag  coefficient  does  not 
necessarily  improve  closure  of  (9)  and  (10)  during  other  periods,  probably  because  of  real  variability 
in  the  effective  bottom  roughness,  which  is  captured  by  the  log-profile  method  but  not  by  a  drag 
law  with  a  fixed  drag  coefficient. 

In  the  discrepancy  in  the  along-isobath  momentum  balance  that  occurs  near  day  10  (Figure 
14),  p0Jo(u-u5)dz  is  much  smaller  than  -Tby.  This  discrepancy  disappears  essentially  completely 
(without  introducing  a  discrepancy  in  the  cross-isobath  momentum  balance)  if  we  set  6  equal  to 
the  height  of  the  top-most  sensor,  rather  than  the  sensor  just  above  the  mixed  layer.  This  fact 
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suggests  that  mixed  layer  thickness  is  a  poor  estimate  of  boundary  layer  thickness  during  this 
particular  event;  i.e.,  that  rzy  just  above  the  mixed  layer  is  not  negligible,  as  required  at  z  =  6  in 
the  mathematical  model  (Section  2.2).  These  ideas  are  amplified  by  an  examination  of  velocity  and 
temperature  profiles.  The  mean  temperature  profile  during  day  10  (Figure  18b)  indicates  a  mixed 
layer  with  a  height  of  approximately  5  m.  Above  the  mixed  layer,  there  is  a  thick  transition  layer 
between  heights  of  roughly  5  and  20  m  above  bottom,  where  the  stratification  is  much  larger  than 
the  stratification  above  20  m.  The  cross-isobath  velocity  is  nonzero  and  sheared  in  the  transition 
layer  as  well  as  the  mixed  layer  (Figure  18a).  It  is  likely  that  turbulent  momentum  transfer  extended 
above  the  mixed  layer  and  well  into  the  transition  layer  during  this  period,  so  that  mixed— layer 
height  is  a  poor  estimate  of  boundary  layer  height.  We  believe  that  the  vertical  structure  in  Figure 
18  is  typical  of  strong  upwelling  events,  but  this  hypothesis  is  difficult  to  test  based  on  the  existing 
observations  because  strong  upwelling  events  are  rare  in  the  STRESS  measurements.  The  vertical 
structure  during  the  more  common  strong  downwelling  events  indicates  a  thick  well-mixed  layer 
overlain  by  a  stratified  interior,  without  a  transition  layer,  so  that  mixed-layer  height  is  probably 
an  accurate  estimate  of  boundary  layer  height  during  downwelling. 

Dever  and  Lentz  (1994)  observed  anomalously  large  horizontal  temperature  gradients  near  day 
40  in  the  SMILE  measurements,  which  had  more  extensive  spatial  coverage  than  the  STRESS 
measurements,  and  they  tentatively  attributed  this  behavior  to  the  presence  of  a  mesoscale  eddy 
on  the  shelf.  It  is  likely  that  during  this  period  the  scales  of  the  density  field  changed  and  that  our 
estimates  of  cross-isobath  and  along-isobath  buoyancy  forces,  which  are  crude  and  involve  implicit 
assumptions  about  spatial  scales,  are  inaccurate.  Failure  to  estimate  buoyancy  forces  accurately 
might  have  caused  the  discrepancies  in  the  momentum  balances  that  occur  near  day  40. 

While  not  definitive,  the  above  discussion  suggests  that  shortcomings  in  our  ability  to  close 
the  momentum  balances  based  on  the  STRESS— 1  measurements  result  from  occasional  problems 
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in  estimation  of  bottom  stress,  boundary  layer  thickness,  and  cross-isobath  and  along-isobath 
buoyancy  forces.  During  STRESS-2,  failures  to  close  the  momentum  balances  are  less  event-like 
(Figures  13  and  15)  and  not  as  clearly  linked  to  specific  explanations,  but  similar  problems  probably 
occurred.  It  is  noteworthy  that  our  results  do  not  suggest  shortcomings  in  the  dynamics  represented 
by  (9)  and  (10),  but  instead  in  our  ability  to  estimate  the  terms  in  these  equations  based  on  the 
existing  measurements. 


5.2  A  balance  for  temperature  anomaly 


The  heat  balance  (9)  is  useful  because  it  permits  us  to  examine  the  relative  importance  of  cross¬ 
isobath  and  along-isobath  advection  in  determining  the  vertically  integrated  variability  of  tem¬ 
perature  within  the  boundary  layer.  A  more  dynamically  relevant  balance,  however,  involves  the 
vertically  integrated  anomaly  in  temperature  within  the  boundary  layer,  relative  to  the  temperature 
in  the  overlying  flow.  As  noted  in  Section  2.3,  temperature  anomaly  is  more  relevant  dynamically 
than  temperature  because  the  anomaly  T  —  T$  appears  in  the  buoyancy  term  in  the  cross— isobath 
momentum  balance  and  because  it  is  related  to  stratification,  which  influences  vertical  mixing. 


To  obtain  a  vertically  integrated  balance  for  temperature  anomaly  that  is  sufficiently  simple 
for  an  observational  test,  let  h  denote  a  fixed  height  large  enough  that  z  =  h  is  always  above  the 
boundary  layer,  so  that  turbulent  fluxes  of  heat  and  momentum  are  negligible  at  z  =  h.  As  in 
Section  2,  we  assume  that  w  dT/dz  udT/dx  and  that  dT/dx  and  dT/dy  are  independent  of  z. 
Equation  (4)  then  yields 


dTh 

dt 


8T  dT 

+  uh—  +  vh—  =  0, 
ox  oy 


(13) 


where  subscript  h  denotes  evaluation  at  z  —  h.  By  subtracting  (13)  from  (4),  integrating  the 
resulting  expression  from  z  =  0  to  z  =  h,  and  applying  the  condition  <j>  =  0  at  both  limits,  one 
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obtains 

i/>-^=-s  sy^-rjy^’ 

wliicb  is  a  vertically  integrated  balance  for  the  anomaly  T  —  2\.  The  term  on  the  left  side  is  the 
time  derivative  of  the  vertically  integrated  temperature  anomaly,  and  the  terms  on  the  right  side 
represent  the  effects  of  cross-isobath  and  along-isobath  advection.  Note  that  in  (13)  we  can  change 
the  order  of  evaluation  at  z  =  h  and  differentiation  with  respect  to  t  because  h  is  fixed.  Similarly, 
in  (14),  use  of  fixed  h  permits  changing  the  order  of  integration  with  respect  to  z  and  differentiation 
with  respect  to  t.  Use  of  the  variable  height  z  =  S(t )  would  have  resulted  in  additional  terms  in 
(14),  which  are  extremely  difficult  to  evaluate  observationally. 

Our  attempts  to  close  (14)  based  on  the  STRESS  measurements,  with  h  fixed  at  the  height 
of  the  top-most  sensor  or  any  other  sensor,  fail  convincingly.  The  magnitudes  of  our  estimates 
of  the  advective  terms  in  (14)  are  comparable  to  that  of  the  time-derivative  term,  but  the  means 
are  different  and  the  variability  is  essentially  uncorrelated.  The  squared  correlation  coefficient  r2 
between  the  time  derivative  term  and  the  sum  of  the  advective  terms  is  0.01  for  both  STRESS-1 
and  STRESS-2,  which  is  well  below  the  critical  value  required  to  reject  the  null  hypothesis  of 
zero  correlation  at  the  95%  confidence  level.  These  results  do  not  change  if  we  neglect  one  of 
the  advective  terms,  and  they  also  do  not  change  if  we  assume  that  the  subinertial  momentum 
balances  are  exactly  correct,  so  that  we  can  use  the  right  sides  of  (9)  and  (10)  (with  S  replaced  by 
h)  to  estimate  one  or  both  of  Jq(u  —  Uh )  dz  and  f{/(v  —  Vh)dz.  Failure  to  close  (14)  is  probably 
a  consequence  of  the  fact  that  the  terms  in  this  balance  involve  only  differences  in  velocity  and 
temperature,  which  are  perhaps  not  accurately  determined  by  the  measurements,  and  the  fact  that 
the  advective  terms  involve  dT/dx  and  dT/dy,  which  are  not  well  resolved  by  the  measurements.  A 
related  issue  is  that  the  assumption  of  vertically  uniform  dT/dy  is  probably  unrealistic,  particularly 
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above  the  bottom  mixed  layer. 


5.3  Other  features  of  the  measurements 


One  of  the  most  striking  features  of  the  observations  is  the  mean  cross-isobath  velocity  profile 
(Figure  10a),  which  indicates  an  unsheared  interior  above  a  sheared  near-bottom  flow  with  a 
thickness  of  approximately  30  m.  To  discuss  this  profile,  note  that  in  (10)  we  can  replace  the 
variable  boundary  layer  thickness  S(t)  with  the  fixed  height  h  introduced  in  Section  5.2,  here  set 
equal  to  the  height  of  the  top-most  sensor  in  the  STRESS  array.  The  mean  of  the  resulting 
expression  is 


rn. _  h2  dT  _ 

Pof  Jo  (u  -  uh)  dz  =  P9-j~q^  ~  Tby 


(15) 


where  an  overbar  denotes  a  mean  value.  The  quantity  on  the  left  side  of  (15)  is  larger  than  the 
stress  on  the  right  side  by  a  factor  of  approximately  eight  for  STRESS-1  and  four  for  STRESS-2 
(Table  1),  so  that  stress  plays  a  minor  role.  Therefore,  if  we  expand  our  viewpoint  from  considering 
only  the  bottom  mixed  layer  (z  <  6)  to  considering  the  bottom  few  tens  of  meters,  we  find  that  the 
dominant  forcing  producing  a  mean  cross-isobath  transport,  relative  to  the  overlying  unsheared 
flow,  is  not  bottom  stress,  as  is  classically  assumed,  but  must  instead  be  some  other  mechanism. 


The  measurements  do  not  indicate  clearly  what  balances  the  left  side  of  (15).  Our  estimate  of  the 
buoyancy  term  on  the  right  side  (obtained  directly  only  during  STRESS-1)  is  smaller  by  a  factor  of 
2.4  than  the  left  side  (Table  1),  so  that  the  balance  does  not  close.  This  observation  is  interesting 
because  the  mean  of  the  corresponding  expression  (10)  closes  quite  well  for  STRESS-1  (Table 
1).  A  significant  nonlinear  effect  does  not  seem  likely  because  of  the  relatively  long  cross-isobath 
and  along-isobath  scales  of  the  topography  (Figure  2)  and  because  crude  estimates  of  vdv/dy 
based  on  differences  between  C3  and  C3'  (possible  only  for  STRESS-1)  suggest  that  this  term 
makes  a  contribution  much  smaller  than  the  left  side  of  (15).  Our  opinion  is  that  our  estimates 
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of  the  buoyancy  term  on  the  right  side  of  (15)  are  not  particularly  accurate,  particularly  above 
the  boundary  layer,  where  our  assumption  that  dT/dy  is  independent  of  z  might  be  particularly 
unrealistic,  and  that  buoyancy  is  in  fact  the  dominant  effect  balancing  the  left  side  of  (15).  This 
idea  is  interesting  because  the  buoyancy  force  must  be  trapped  in  the  bottom  30  m  in  order  to 
explain  the  unsheared  cross-isobath  velocity  above  30  m. 

An  additional  striking  feature  of  the  observations  is  the  Ekman-like  structure  indicated  by 
the  most  energetic  EOF  of  velocity  (Figure  10),  which  represents  all  but  a  small  fraction  of  the 
variance  in  the  subinertial  velocity  field.  It  is  natural  to  ask  whether  this  EOF  representation  is  more 
consistent  with  the  vertically  integrated  momentum  balances  than  are  the  original  measurements. 
We  have  found  that  the  time  series  described  by  the  most  energetic  EOF  mode  for  velocity  and 
most  energetic  EOF  mode  for  temperature  are  in  substantially  worse  agreement  with  (9)  and  (10) 
than  are  the  original  time  series.  Thus  the  velocity  and  temperature  profiles  represented  by  the 
most  energetic  EOF  modes  do  not  reflect  the  dynamics  represented  by  the  vertically  integrated 
momentum  balances. 

A  brief  consideration  of  temporal  acceleration,  which  was  neglected  in  the  vertically  inte¬ 
grated  momentum  balances,  is  convenient  at  this  point.  If  one  includes  temporal  acceleration 
in  (1)  and  (2)  and  uses  Leibniz’s  rule  to  interchange  the  order  of  integration  with  respect  to 
z  and  differentiation  with  respect  to  t  in  deriving  (7)  and  (8)  from  (5)  and  (6),  the  result  is 
addition  of  d/dt  /05  p0(u  -  us)dz  +  poS(d8/dt)(du/dz)s  to  the  left  side  of  (9)  and  addition  of 
d/dt  /05  po(v  —  vg)dz  +  po8(dS/dt)(dv/dz)s  to  the  left  side  of  (10).  We  have  confidence  in  our 
ability  to  estimate  d/dt  Jq  ( u  —  ug)dz  and  d/dt  Jq(v  —  v$)dz  with  reasonable  accuracy  based  on 
the  STRESS  measurements,  but  we  have  no  confidence  in  estimates  of  the  terms  involving  d6/dt , 
primarily  because  these  terms  are  highly  nonlinear  and  because  of  the  coarseness  of  the  STRESS 
array  at  heights  above  6  m  (Figure  3),  which  makes  evaluation  of  velocity  gradients  and  dS/dt 
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problematic.  Consequently,  a  quantitative  assessment  of  the  importance  of  temporal  acceleration  is 
not  possible.  A  qualitative  assessment  can  be  obtained,  however,  by  considering  only  the  quantities 
d/dt  fo  p0(u  -  u$)  dz  and  d/dt  /05  pQ(v  -  vs)dz.  The  means  of  these  quantities  are  two  orders  of 
magnitude  smaller  than  the  means  of  the  leading  terms  in  (9)  and  (10),  and  their  standard  deviation 
are  smaller  by  at  least  a  factor  of  five  than  the  standard  deviations  of  the  leading  terms  (Table  1). 
These  results  indicate  that  temporal  acceleration  has  only  a  minor  effect  on  vertically  integrated 
momentum  balances  for  the  subinertial  flow.  This  conclusion  is  not  surprising  for  the  cross-isobath 
momentum  balance,  since  du/dt  <C  f  v  for  sub-inertial  flows,  but  it  is  at  first  surprising  for  the 
along- isobath  momentum  balance,  since  dv/dt  and  fu  are  comparable.  The  reason  that  temporal 
acceleration  is  not  important  in  the  vertically  integrated  along-isobath  momentum  balance  is  that 
( u  —  u$)  and  (v  -  vg)  are  comparable,  even  though  u  <  v,  so  that  d/dt(v  -  vs)  <  f(u  -  us)  for 
sub-inertial  fluctuations. 

6  Conclusion 

We  have  presented  an  observational  test  of  a  set  of  vertically  integrated,  subinertial  balances  for 
momentum  and  heat  in  the  bottom  boundary  layer,  based  on  time-series  measurements  obtained 
on  the  northern  California  shelf  during  each  of  two  winter  seasons.  As  has  recently  been  predicted 
theoretically,  the  measurements  indicate  clearly  that  the  cross-isobath  momentum  equation  does 
not  reduce  simply  to  a  classical  Ekman  balance,  but  instead  includes  a  dynamically  significant 
buoyancy  force  resulting  from  the  distortion  of  isopycnal  surfaces  within  the  boundary  layer.  Also 
as  predicted  theoretically,  the  cross-isobath  buoyancy  force  is  dominant  when  the  boundary  layer 
is  thick,  which  typically  occurs  during  downwelling-favorable  flows.  Subinertial  fluctuations  in 
the  along-isobath  momentum  balance  are  crudely  consistent  with  an  Ekman  balance.  The  mean 
along-isobath  momentum  balance  is  significantly  affected  by  a  buoyancy  force  produced  by  an 
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along-isobath  temperature  gradient,  and  along-isobath  advection  is  as  important  as  cross-isobath 
advection  in  the  heat  balance. 

Two  important  conclusions  result  from  our  work.  The  first  is  that  a  buoyancy  force  produced 
by  distortion  of  the  isopycnal  surfaces  within  the  boundary  layer  is  a  significant  term  in  the  cross- 
isobath  momentum  balance,  and  that  the  buoyancy  force  behaves  in  approximately  the  manner 
predicted  by  the  recent  theoretical  models.  The  second  conclusion  is  that  along-isobath  variability 
of  temperature  has  an  important  effect  on  the  bottom  boundary  layer.  The  along-isobath  temper¬ 
ature  gradient  influences  both  the  heat  balance,  where  it  contributes  along-isobath  advection,  and 
the  mean  momentum  balance,  where  it  produces  an  along-isobath  buoyancy  force  that  modifies  the 
cross-isobath  velocity.  Thus  along-isobath  variability  of  temperature  influences  both  the  evolution 
of  the  temperature  field  and  the  structure  of  the  velocity  field.  These  effects  have  been  neglected 
in  the  theoretical  work  on  the  bottom  boundary  layer  that  was  reviewed  by  Garrett  et  al.  (1993) 
and  in  most  numerical  simulations  of  the  entire  shelf  (e.g.,  Allen  et  al.,  1995;  Allen  and  Newberger, 
1996),  although  they  have  been  included  in  the  theoretical  model  developed  recently  by  Chapman 
and  Lentz  (in  press). 

An  important  observation  also  results  from  our  work.  The  measurements  indicate  a  mean  cross¬ 
isobath  velocity  that  is  offshore  below  20  m  above  bottom,  onshore  above  20  m,  and  uniform  above 
30  m.  The  transport  produced  by  this  flow,  relative  to  the  unsheared  flow  above  30  m,  is  larger  by 
a  factor  of  four  to  eight  than  the  Ekman  transport  associated  with  the  mean  along-isobath  bottom 
stress.  Thus  the  dominant  mechanism  producing  a  mean  cross-isobath  transport,  relative  to  the 
unsheared  flow  above  30  m,  is  not  bottom  stress,  as  is  classically  assumed,  but  must  be  some  other 
mechanism.  A  likely  mechanism  is  a  vertically  nonuniform  along-isobath  buoyancy  force,  which  is 
poorly  resolved  by  the  measurements. 
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Table  1:  Statistics  of  the  Terms  in  the  Vertically  Integrated  Balances  (SI  Units) 
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mean  std  dev 

mean  std  dev 
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'cT 

1 

O", 

+0.0279  0.0338 

+0.0315  0.0518 

+0.0100  0.0198 

+0.0149  0.0249 

Pgl(P/2)(dT/dx)  +  a  }‘(T  -  T, )  dz] 

+0.0113  0.0303 

+0.0208  0.0549 

Po(d/ dt)  Jq(u  -  us)dz 

-0.0001  0.0051 

+0.0001  0.0067 

-0.0114  0.0257  -0.0209  0. 

-0.0064  0.0305  -0.0124  0. 

-0.0039  0.0057  -0.0059 

-0.0001  0.0060  +0.0001  0. 

+0.1128  0.3417  +0.0676 

+0.0641  0.2734  +0.1063 

+0.1345  0.2704  +0.1580 


-0.0498 

-0.0124 

-0.0198 
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Figure  Captions 

Figure  1.  Configuration  addressed  in  the  theoretical  development.  A  non-conducting  sea  floor 
is  inclined  at  a  locally  uniform  angle  a  with  respect  to  horizontal.  A  turbulent  boundary  layer 
adjacent  to  the  sea  floor  is  overlain  by  a  stably  stratified  interior  in  which  turbulent  fluxes  of 
momentum  and  heat  are  negligible. 

Figure  2.  STRESS  study  area  on  the  northern  California  shelf.  C3  is  the  central  site  and  C2,  C3' 
and  C4  are  secondary  sites.  C3  and  C3'  are  at  water  depths  of  approximately  90  m  and  C4  is  at  a 
water  depth  of  approximately  130  m.  The  coordinate  system  is  oriented  so  that  x  is  cross-isobath 
and  y  is  along-isobath  at  C3. 

Figure  3.  STRESS  instrumentation.  During  STRESS-1,  the  heights  of  the  BASS  sensors  were 
0.21,  0.51,  1.10,  1.96,  2.56  and  4.96  m;  the  heights  of  the  temperature  sensors  were  0.37,  0.62,  1.10, 
2.01,  2.51,  3.63,  5.00  and  5.84  m;  and  the  heights  of  the  VACMs  were  6,  12, 18,  24  and  30  m  at  C3 
and  6  and  18  m  at  C3'.  During  STRESS-2,  the  heights  of  the  BASS  sensors  at  C3  and  C4  were 
0.39,  0.74,  1.34,  1.94,  2.54  and  4.94  m;  the  heights  of  the  temperature  sensors  were  0.19,  1.08,  1.75, 
2.62,  3.76,  4.48,  5.14  and  5.83  m  at  C3  and  0.15,  2.93,  3.77,  4.53,  5.18,  5.96,  13,  23,  27,  41  and  56 
m  at  C4;  the  heights  of  the  temperature-conductivity  sensors  were  16,  22,  25  and  28  m  at  C3  and 
7  m  at  C4;  and  the  heights  of  the  VACMs  were  7,  10,  13,  19  and  31  m  at  C3  and  10,  19  and  71 
m  at  C4.  During  STRESS-2,  the  third  and  sixth  BASS  above  bottom  at  C3,  the  first  and  third 
BASS  above  bottom  at  C4,  the  second  VACM  velocity  sensor  above  bottom  at  C3,  and  the  lowest 
VACM  conductivity  sensor  at  C3  malfunctioned,  as  noted  in  the  text.  We  do  not  use  the  highest 
STRESS-2  VACM  record  at  C4,  which,  at  71  m  above  bottom,  was  far  above  the  bottom  mixed 
layer.  Sensors  above  35  m  are  not  shown  in  the  figure. 

Figure  4.  Temperature-density  relationships  based  on  hourly  data  at  C3.  The  STRESS-1  mea- 
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s u.rement s  were  obtained  from  temperature— conductivity  sensors  at  heights  of  6, 18  and  30  m  above 
bottom.  The  STRESS-2  measurements  were  obtained  from  temperature-conductivity  sensors  at 
heights  of  10,  13,  16,  19,  22,  25,  28  and  31  m.  The  solid  lines  represent  least-squares  fits  to  the 
temperature-density  relationship,  which  correspond  to  j3  =  0.35  kg  m-3  deg-1  C  for  STRESS-1 
and  j3  =  0.29  kg  m-3  deg-1  C  for  STRESS-2. 

Figure  5.  Relationship  between  current  speed,  obtained  from  velocity  averaged  over  the  bottom 
three  BASS  sensors,  and  log-profile  estimates  of  shear  velocity  for  STRESS-1  (left  panel)  and 
STRESS-2  (right  panel)  based  on  the  hourly  data.  The  solid  lines  correspond  to  a  quadratic  drag 
law  with  a  bottom  roughness  parameter  of  zq  =  0.01  m. 

Figure  6.  Rotary  velocity  spectra  based  on  all  functioning  current  sensors  from  STRESS-1  (top 
panel)  and  STRESS-2  (bottom  panel).  Positive  and  negative  frequencies  correspond  to  counter¬ 
clockwise  and  clockwise  polarization,  respectively.  The  inertial  frequency  is  approximately  -1.25 
cpd.  We  estimated  spectra  by  detrending  and  demeaning  the  data,  applying  a  Hanning  window, 
computing  the  discrete  Fourier  transform  of  the  data  in  the  complex  form  u  +  iv,  squaring  the 

absolute  value  of  the  Fourier  coefficients,  and  averaging  over  17  adjacent  frequencies  to  obtain 

spectral  estimates  with  34  degrees  of  freedom.  We  obtained  spectra  for  STRESS-1  by  concatenat¬ 
ing  the  time  series  for  STRESS-la  and  STRESS-lb.  The  more  energetic  spectra  correspond  to 
measurements  higher  in  the  water  column. 

Figure  7.  Rotary  spectra  of  bottom  stress  during  STRESS-1  and  STRESS-2,  computed  as  in 
Figure  6. 

Figure  8.  Time  series  of  filtered  velocity  and  temperature  during  STRESS-1. 

Figure  9.  Time  series  of  filtered  velocity  and  temperature  during  STRESS-2. 
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Figure  10.  Mean  velocity  and  temperature  together  with  EOF  representation  of  filtered  velocity 
and  temperature  fluctuations  during  STRESS-1  and  SMILE  (circles)  and  STRESS-2  (pluses). 
Measurements  above  40  m  are  from  the  SMILE  array. 

Figure  11.  Estimates  of  gradient  Richardson  number  Ri  at  four  heights  based  on  hourly  STRESS- 
1  measurements.  The  gradients  are  calculated  based  on  differences  between  measurements  at  two 
different  heights.  The  heights  listed  for  Ri  are  the  means  of  the  heights  of  each  pair  of  measurements. 
The  pairs  of  heights  are  18  and  12  m  for  the  top  panel,  12  and  6  m  for  the  second  panel,  4.96  and 
2.56  m.for  the  third  panel,  and  2.56  and  0.21  m  for  the  fourth  panel.  Note  that  the  vertical  spacing 
between  the  pair  of  sensors  is  approximately  2.5  m  for  the  bottom  two  panels  and  6  m  for  the  top 
two  panels.  The  solid  lines  correspond  to  Ri  =1/4. 

Figure  12.  Terms  in  the  cross-isobath  momentum  balance  (9),  estimated  from  the  STRESS-1 
measurements.  For  the  top  panel,  r2  -  0.17,  and  for  the  center  panel,  r2  =  0.55.  For  convenience, 
the  bottom  panel  shows  the  thickness  of  the  bottom  mixed  layer. 

Figure  13.  Terms  in  the  cross-isobath  momentum  balance  (9),  estimated  from  the  STRESS-2 
measurements.  For  the  top  panel,  r2  =  0.15,  and  for  the  center  panel,  r2  =  0.54.  For  convenience, 
the  bottom  panel  shows  the  thickness  of  the  bottom  mixed  layer. 

Figure  14.  Terms  in  the  along-isobath  momentum  balance  (10),  estimated  from  the  STRESS-1 
measurements.  For  the  top  panel,  r2  =  0.33,  and  for  the  bottom  panel,  r2  =  0.33. 

Figure  15.  Terms  in  the  along-isobath  momentum  balance  (10),  estimated  from  the  STRESS-2 
measurements.  For  the  top  panel,  r2  =  0.26,  and  for  the  bottom  panel,  r2  —  0.34. 

Figure  16.  Terms  in  the  integrated  heat  balance  (11),  estimated  from  the  STRESS-1  measurements. 
For  the  top  panel,  r2  =  0.11,  and  for  the  bottom  panel,  r2  =  0.35. 
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Figure  17.  Terms  in  the  integrated  heat  balance  (11),  estimated  from  the  STRESS-2  measurements. 
For  the  top  panel,  r 2  =  0.27,  and  for  the  bottom  panel,  r2  =  0.50. 

Figure  18.  Mean  velocity  and  temperature  profiles  for  upwelling  event  near  day  10  during  STRESS- 1 
(January  10,  1989). 
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